Draft version February 5, 2008 

Preprint typeset using I4TgX style emulateapj v. 6/22/04 



MAGNETIC HELICITY DENSITY AND ITS FLUX IN WEAKLY INHOMOGENEOUS TURBULENCE 

Kandaswamy Subramanian 1 and Axel Brandenburg 2 

Draft version February 5, 2008 

ABSTRACT 

A gauge invariant and hence physically meaningful definition of magnetic helicity density for random fields 
is proposed, using the Gauss linking formula, as the density of correlated field line linkages. This definition 
is applied to the random small scale field in weakly inhomogeneous turbulence, whose correlation length is 
small compared with the scale on which the turbulence varies. For inhomogeneous systems, with or without 
boundaries, our technique then allows one to study the local magnetic helicity density evolution in a gauge 
independent fashion, which was not possible earlier. This evolution equation is governed by local sources 
(owing to the mean field) and by the divergence of a magnetic helicity flux density. The role of magnetic 
helicity fluxes in alleviating catastrophic quenching of mean field dynamos is discussed. 

Subject headings: MHD — turbulence — Sun: magnetic fields 



1. INTRODUCTION 

Large scale magnetic fields produced by dynamo action 
tend to have some degree of magnetic helicity. A simple ex- 
ample is the interlocking of poloidal and toroidal fields in one 
hemisphere of a star, seen from stellar dynamo models. How- 
ever, in the case of the Sun there is explicit evidence of mag- 
netic helicity being present in or coming from active regions 
(Seehafer 1990; Pevtsov et al. 1995; Bao et al. 1999; Berger 
& Ruzmaikin 2000), coronal mass ejections (Demoulin et al. 
2002), and the solar wind (Matthaeus et al. 1982; Lynch et 
al. 2005). While the investigation of magnetic helicity in the 
Sun and in the solar wind is interesting in its own right, there 
is now also a direct connection with dynamo theory with the 
realization that large scale dynamos must transport and shed 
small scale magnetic helicity in order to operate on a dynami- 
cal time scale rather than the much longer resistive time scale 
(see the review of Brandenburg & Subramanian 2005a for ref- 
erences). A major difficulty with this picture is the absence 
of a meaningful definition for magnetic helicity density, even 
for small scale fields. Magnetic helicity is a volume integral, 
usually defined as Hm — f A ■ B dV, where A is the vector 
potential and B = V x A is the magnetic field. However, 
under a gauge transformation A' = A + VA, which leaves 
B invariant, one has H' = H + J KB ■ dS. So H is only 
gauge invariant if the B field has no component normal to the 
boundary or if it vanishes sufficiently rapidly at the bound- 
ary of the integration volume. In most practical contexts, like 
the Sun or galaxies, however, the field does not vanish on the 
boundaries. A possible remedy might be to consider instead 
the gauge-invariant relative magnetic helicity, defined by sub- 
tracting the helicity of a reference vacuum field in the same 
gauge (Berger & Field 1984, Finn & Antonsen 1985). But 
the flux of relative helicity is cumbersome to work with for 
arbitrarily shaped boundaries. Also the concept of a density 
of relative helicity is not meaningful, since it is defined only 
as a volume integral. Indeed, there is simply no way that the 
quantity A ■ B itself can be gauge-invariant. 

On an earlier occasion, Subramanian & Brandenburg 
(2004, hereafter SB04) considered the evolution of the cur- 
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rent helicity density, Hq = J ■ B, where J = V x B/ fir, 
is the current density and /in. is the vacuum permeability (we 
set no = 1 in what follows). Note that Hq, as well as its flux, 
are locally well defined, explicitly gauge invariant and obser- 
vationally measurable. Furthermore, from a closure model, 
Pouquet et al. (1976) show that the a effect needed for large- 
scale dynamos has a nonlinear addition due to the small scale 
contribution to He- The build up of this small-scale current 
helicity then goes to cancel the kinetic part of the a effect and 
causes catastrophic quenching of the dynamo, unless one can 
have a helicity flux out of the system. This formed the motiva- 
tion for the work of SB04. The major disadvantage in work- 
ing with Hq however is that one loses the conceptually simple 
form of the magnetic helicity conservation law. We propose 
here instead an alternative means to define magnetic helicity 
density for the random small scale field, using the more basic 
Gauss linking formula for helicity, which can be directly ap- 
plied to discuss magnetic helicity density and its flux even in 
inhomogeneous systems with boundaries. The technique ap- 
plied in calculating the magnetic helicity evolution is however 
very similar to that employed in SB04. 

In the following we define random small scale quantities as 
departures from the corresponding mean field quantity, e.g., 
b = B — B for the magnetic field, j = J — J for the cur- 
rent density, and u = U — U for the velocity. Throughout 
this Letter we adopt ensemble averages that, in practice, are 
commonly approximated as spatial averages over one or two 
coordinate directions (see, e.g., Brandenburg & Subramanian 
2005a). However, the approach developed below applies also 
to the case without a mean magnetic field. Therefore, specific 
applications to the mean field dynamo (MFD) will be post- 
poned until the end of this Letter. 



2. MAGNETIC HELICITY DENSITY 

Given the random small scale magnetic field b(x,t) one 
can also define the magnetic helicity directly in terms of the 
field, as the linkage of its flux, using Gauss's linking formula 
(Berger & Field 1 984, Moffatt 1 969) 
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where both integrations extend over the full volume. Suppose 
we define an auxiliary field 
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then this field satisfies V x a = b, and V ■ a = 0, and one 
can write Iiq = jab d 3 x. This is the origin of the text- 
book definition of magnetic helicity in what is known as the 
Coulomb gauge for the vector potential. Provided the field 
is closed over the integration volume, this definition can be 
applied in any other gauge. Note however that for an open 
system with boundaries, it is not useful to go to the definition 
involving the vector potential, which is now of course gauge 
dependent. We therefore take the point of view here that the 
magnetic helicity density Hq defined by equation Q is the 
more basic definition of the topological property determining 
the links associated with magnetic fields, and not the defini- 
tion in terms of the vector potential. We will see below that 
this then also allows us to define naturally a gauge invariant 
magnetic helicity density for random fields, as long as the cor- 
relation scale of the field is much smaller than the size of the 
system, as the density of correlated links of the field. 

For th is consider b to be a random field with a correlation 
function bi(x,t)bj(y,t) — Mij(r,R). Here we have defined 
the difference r = x — y and the mean R = (x + y)/2, 
keeping in mind that, for weakly inhomogeneous turbulence, 
the two-point correlation M,; 3 (r, R) and in fact all two-point 
correlations below, vary rapidly with r but vary slowly with 
R. Taking the ensemble average of Hq, we have 
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Next, we suppose that the correlation scale I of the random 
small scale field b is much smaller than the system scale i?g, 
i.e. we suppose that there exists an intermediate scale L such 
that I < L < R s with Mij(r,R) -» as \r\ -> L > 
I. Then one can do the r integral even by restricting oneself 
to the intermediate scale L and still capture all the dominant 
contributions to the integral. This then motivates us to define 
the magnetic helicity density h of the random small scale field 
as 7T G = / d 3 R h(R). Here, 
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where we can formally let L — > oo. The above expression for 
h(R) in equation (0} is our proposal for the helicity density of 
the random small scale field b. Evidently, h(R) is explicitly 
gauge invariant. A qualitative description would be to say that 
the magnetic helicity density of a random small scale field is 
the density of correlated links of the field. We can now de- 
rive the evolution equation for h(R) and also meaningfully 
(in a gauge invariant manner) talk about its flux. Note that 
this has not been possible before, although many papers (e.g., 
Blackman & Field 2000; Kleeorin et al. 2000; Vishniac & Cho 
2001) have appealed to the notion of a magnetic helicity flux 
density in some qualitative fashion. We will see that the mag- 
netic helicity evolution equation that we derive reproduces the 
known evolution equation for homogeneous turbulence and 
generalizes it to the inhomogeneous case by introducing pos- 
sible fluxes of helicity. 

3. MAGNETIC HELICITY DENSITY EVOLUTION 



It is much simpler to work out the evolution equation for 
h(R) by first going to Fourier space, using the two-scale ap- 
proach of Roberts & Soward (1975), where all two-point cor- 
relations are assumed to vary rapidly with r and slowly with 
R. Consider the equal time, ensemble average of the prod- 
uct f(xi)g(x2)- The common dependence of / and g on 
t is assumed and will not explicitly be stated. Let /(fei) 
and g(k2) be the Fourier transforms of / and g, respec- 
tively. We can express this correlation as f(xi)g(x2) — 
J $(/, g, k, R) e ik r d 3 k, with 



$(f,g,k,R) = J f{k + \K)g(-k + \K) 
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d 6 K. 
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Here, k — h(ki — £2) andK = ki+k2- We define in Fourier 
space the correlation and cross correlation tensors of the u 
and b fields; Vij(k,R) = 5>(uj, uj, k, R), rriij(k,R) = 

$(bi,bj,k,R), and Xjk{k,R) = &(y,j,bk,k,R). In 
MFD theory, the turbulent electromotive force (EMF) is 
given by £ = u x b, whose components are £i(R) = 
djk J Xjk{k, R) d 3 fc. Furthermore, in Fourier space we have 
for the magnetic helicity density, 



h(R) = I I e ijk h{k + \K%{-k + \K) 
x{ik k /k 2 )e iKR d 3 kd 3 K. 
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We should remark, that for an inhomogeneous system, the 
Coulomb gauge magnetic helicity density, say h = a ■ b, 
would have (k k + Kk/2)/(k + K/2) 2 replacing k^/k 2 in 
the Fourier space expression of equation 0. The two expres- 
sions are identical for the homogeneous case, and even in the 
weakly inhomogeneous case up to first order terms in K/k, 
but not in general. So, h(R) ^ h(R) in general. 

In order to compute the magnetic helicity evolution, we use 
the induction equation for b in Fourier space, dbi(k)/dt = 
—Ci pq \k p e q . Here e q is the Fourier transform of the small scale 
electric field e q , which is given by (e.g., Moffatt 1978) 



e = —u x B — U x b — u x b + u x b + r\j . (7) 

Substituting this in the time derivative of equation (|6}, we get 
after some straightforward algebra, 
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We denote the integrals over the three terms in curly brackets 
above as Ai, A%, and A3, respectively. From the definition of 
$ in Eq. (13, the first term is simply A\ — — 2e ■ b, or 



A± — 2b- (u x B) — 2r)j ■ b = -2£ • B - 2rjj ■ b, (9) 



where £ = u x b is the turbulent EMF. Note that for homo- 
geneous turbulence only the term Ai survives. This is because 
the other two terms which involve a large K t in the integrand, 
will introduce a large scale Ri derivative when evaluating the 
integral over K , and this vanishes in the homogeneous case. 
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So, for the homogeneous case we recover a local general- 
ization (without volume integration) of the magnetic helicity 
conservation equation (Brandenburg & Subramanian 2005a) 
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where the flux is J-s = J- s 
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dh/dt = -25 B-2r]j b (from A 1 term only). (10) 

In the inhomogeneous case, since the terms A 2 and A 3 are 
scalars that depend on Ki in the integrand, and hence a large 
scale Ri derivative, they will contribute purely to the flux of 
helicity. The only term that involves volume generation of the 
helicity density is the A\ term, which we see involves corre- 
lations no higher than the two-point one. This is in contrast to 
the current helicity evolution, which involved undetermined 
triple correlations in their volume generation; see SB04. 

Let us now evaluate the helicity fluxes given by A 2 and A3 . 
This involves straightforward but tedious algebra. We also 
work out the flux to the lowest order in the R derivative. There 
are again three main types of contributions due to different 
parts of the electric field e. First there is a contribution propor- 
tional to B due to that part e = -mxB+... . In Fourier space 

this gives e q (k) = e q i m J u m {k — k')Bi(k')d 3 k' . We substi- 
tute this into the expressions for A 2 and A3, change the vari- 
ables to K' = K k' , use the definition for Xij an d evaluate 
the integrations over K' and k' retaining only terms to lowest 

order in the R derivatives. We then get A 2 



A3 = —VjJ 7 ^, where the mean field dependent fluxes T 
and are given by 
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Note that only depends on the antisymmetric part of the 

— vc 

cross correlation Xmq, whereas T i is sensitive to the sym- 
metric part as well. Now consider the contribution propor- 
tional to the mean velocity from the part of the electric field 
e = — U x b + ... . The evaluation of this follows the same 
steps as in evaluating (II 1> except that one can map u m — > bi 

and Bi — » U m . This gives A 2 + A3 = —VjJ-j , where the 
flux due to bulk motions is given by 

^ Ulk = t qlm U rn {R) J (2ik q mii~ik i m lq )k^ 2 A i k. (12) 

Indeed, if the magnetic correlations were isotropic then it is 

easy to simplify this further and one gets T i U = hUi, ex- 
actly as one should for an advective flux! 

The contribution to the fluxes from e = u x b (see equa- 
tion |7}, introduces triple correlations in the flux, which then 
need a closure theory to evaluate. We denote this flux term as 

•7-"/ lp e . However, since this triple correlation comes only in 
the flux, and not the volume generation term, it is likely that 
its value can be constrained by a conservation law. This will 
be examined in more detail in the future. Note a lso that the 
contribution to the helicity flux from e = u x b + . . is zero 
and that the resistive contribution from e = rjj + ... is likely 
to be negligible compared to the terms that we retain. 

Putting all our results together we can write for the evolu- 
tion of the magnetic helicity density 



dh/dt + V • T = -2£ ■ B - 2r)j ■ b, 
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should emphasize that equation Jl 3i is a local magnetic he- 
licity conservation law. If one were unable to define a gauge- 
invariant magnetic helicity density, one would only have an 
integral (global) conservation law, as in previous studies. Fur- 
ther simplification of the helicity fluxes for use in say MFD 
models, requires the evaluation of the turbulent EMF tensor 
Xij, which can be done only under a closure scheme, and will 

— vc 

be presented elsewhere. It turns out that T \ is a general- 
ization of the magnetic helicity flux obtained by Vishniac & 
Cho (2001), which is particularly important in the presence 
of strong shear (Brandenburg & Subramanian 2005b). The 

A 

fluxes used by Kleeorin et al. (2000) can arise from both T i 

and the contributions of the antisymmetric parts of the corre- 
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lations to T \ . The magnetic helicity fluxes also vanish if the 

turbulence is homogeneous. For isotropic but inhomogeneous 
— vc 

turbulence, T i also depends purely on the antisymmetric 

A 

part of Xij, just like T \ . Furthermore, the VC and A-fluxes 
are proportional to two-point correlations (Xij) an d tne mean 
field B, see equation il Ik this is unlike the two-dimensional 
case (Silvers 2006). 

Since the small scale magnetic helicity opposes the kinetic 
part of the a effect (Pouquet et al. 1976), its loss through cor- 
responding magnetic helicity fluxes can alleviate this quench- 
ing effect (Blackman & Field 2000; Kleeorin et al. 2000; 
Vishniac & Cho 2001). Note that the a effect quantifies 
the contribution of £ that is aligned with the mean field, i.e. 
S = aB + ... for the simplest case of a scalar a effect. For 
a closed system equation dlOi applies and, in the stationary 
limit, this predicts £ ■ B = —rjj ■ b which tends to zero as 
77 — > for any reasonable spectrum of current helicity. This 
leads to a catastrophic quenching of the turbulent EMF paral- 
lel to B. In the presence of helicity fluxes, however, we have 
S B = — • J- — rjj ■ b in the stationary limit, and so 
£ ■ B need not be catastrophically quenched. The turbulent 
magnetic helicity fluxes worked out here are therefore crucial 
for the efficient working of the mean field dynamo. Numer- 
ical work in determining the a effect did show a 30-fold in- 
crease in simulations that allowed helicity fluxes to develop 
(Brandenburg & Sandin 2004). However, a more convincing 
demonstration of the importance of helicity fluxes comes from 
a dynamo simulation in the presence of shear showing that 
only with open boundaries can a significant large scale field 
of equipartition field strength develop (Brandenburg 2005). 

4. CONCLUSIONS 

We have proposed here a local gauge invariant definition of 
magnetic helicity density for random fields in weakly inho- 
mogeneous systems, which can also have boundaries. This is 
particularly useful in the context of MFDs since one can then 
meaningfully discuss magnetic helicity fluxes and the local 
effect of Lorentz forces. We have derived an evolution equa- 
tion for the local magnetic helicity density and showed that 
they naturally involve helicity fluxes, which may alleviate the 
problems associated with MFDs (Shukurov et al. 2006). Our 
work therefore lays the conceptual foundation for the many 
discussions of the effects of helicity fluxes already existing in 
the literature and for future explorations. 

Future applications might include the use of equation © to- 
gether with an assumption of isotropy to estimate the spatial 
variation of magnetic helicity density by measuring at least 
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one of the off-diagonal components of Mij(r, R). This type 
of approach has been adopted by Matthaeus et al. (1982) to 
determine the magnetic helicity in the solar wind by measur- 
ing just a time series of an off-diagonal component of M y (r ) 
under the Taylor hypothesis. However, no dependence on the 
large scale coordinate R has been determined. In principle, 
similar ideas could be applied to determine h(R) on the solar 
surface without necessarily having access to the dependence 
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